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Abstract 

The measurement of the jet energy resolution is presented using data recorded with the atlas 
detector in proton-proton collisions at = 7 TeV. The sample corresponds to an integrated luminosity 
of 35 pb ^ Jets are reconstructed from energy deposits measured by the calorimeters and calibrated 
using different jet calibration schemes. The jet energy resolution is measured with two different in 
situ methods which are found to be in agreement within uncertainties. The total uncertainties on these 
measurements range from 20% to 10% for jets within | v| < 2.8 and with transverse momenta increasing 
from 30 GeV to 500 GeV. Overall, the Monte Carlo simulation of the jet energy resolution agrees with 
the data within 10%. 





EPJ manuscript No. 

(will be inserted by the editor) 



Jet energy resolution in proton-proton collisions at v^ = 7 TeV 
recorded in 2010 with the ATLAS detector 

The ATLAS Collaboration 

Address(es) of author(s) should be given 

Received: date / Revised version: date 

Abstract. Tlie measurement of the jet energy resolution is presented using data recorded with the ATLAS detector 
in proton-proton collisions at y/s = 7 TeV. The sample corresponds to an integrated luminosity of 35 pb^'. Jets are 
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1 Introduction 



Precise knowledge of the jet energy resolution is of key im- 
portance for the measurement of the cross-sections of inclusive 
jets, dijets, multijets or vector bosons accompanied by jets [1- 
4], top-quark cross-sections and mass measurements [5], and 
searches involving resonances decaying to jets [6,7]. The jet 
energy resolution also has a direct impact on the determina- 
tion of the missing transverse energy, which plays an impor- 
tant role in many searches for new physics with jets in the fi- 
nal state [8,9]. This article presents the determination with the 
ATLAS detector [10, 1 1] of the jet energy resolution in proton- 
proton collisions at a centre-of-mass energy of y/s — 7 TeV. 
The data sample was collected during 2010 and corresponds to 
35 pb^ ' of integrated luminosity delivered by the Large Hadron 
Collider (LHC) [12] at CERN. 

The jet energy resolution is determined by exploiting the 
transverse momentum balance in events containing jets with 
large transverse momenta (pj). This article is structured as fol- 
lows: Section 2 describes the ATLAS detector. Sections 3, 4 
and 5 respectively introduce the Monte Carlo simulation, the 
event and jet selection criteria, and the jet calibration methods. 
The two techniques to estimate the jet energy resolution from 
calorimeter observables, the dijet balance method [13] and the 
bisector method [14], are discussed respectively in Sections 6 
and 7. These methods rely on somewhat different assumptions, 
which can be validated in data and are sensitive to different 
sources of systematic uncertainty. As such, the use of these two 
independent in situ measurements of the jet energy resolution 
is important to validate the Monte Carlo simulation. Section 8 
presents the results obtained for data and simulation for the 
default jet energy calibration scheme implemented in ATLAS. 
Section 9 compares results of the Monte Carlo simulation in 
situ methods to the resolutions obtained by comparing the jet 
energy at calorimeter and particle level. This comparison will 
be referred to as a closure test. Sources of systematic uncer- 
tainty on the jet energy resolution estimated using the available 
Monte Carlo simulations and collision data are discussed in 
Section 10. The results for other jet energy calibration schemes 
are discussed in Sections 1 1 and 12, and the conclusions can be 
found in Section 13. 



2 The ATLAS detector 

The ATLAS detector is a multi-purpose detector designed to 
observe particles produced in high energy proton-proton col- 
lisions. A detailed description can be found in Refs. [10, 11]. 
The Inner (tracking) Detector has complete azimuthal cover- 
age and spans the pseudorapidity region jTj | < 2.5 ' . The Inner 



^ The ATLAS reference system is a Cartesian right-handed coor- 
dinate system, with the nominal collision point at the origin. The 
anti-clockwise beam direction defines the positive z-axis, with the 
jj-axis pointing to the centre of the LHC ring. The angle (/> defines 
the direction in the plane transverse to the beam {x,y). The pseudora- 
pidity is given by T] = — Intan-j, where the polar angle is taken 
with respect to the positive z direction. The rapidity is defined as 
y = 0.5 X ln[{E + p,) / {E — p~)], where E denotes the energy and is 
the component of the momentum along the z-axis. 



Detector consists of layers of silicon pixel, silicon microstrip 
and transition radiation tracking detectors. These sub-detectors 
are surrounded by a superconducting solenoid that produces a 
uniform 2 T axial magnetic field. 

The calorimeter system is composed of several sub- 
detectors. A high-granularity liquid-argon (LAr) electromag- 
netic sampling calorimeter covers the \rj\ < 3.2 range, and it 
is split into a barrel < 1.475) and two end-caps (1.375 < 
\ri \ < 3.2). Lead absorber plates are used over its full cover- 
age. The hadronic calorimetry in the barrel is provided by a 
sampling calorimeter using steel as the absorber material and 
scintillating tiles as active material in the range |t7| < 1.7. This 
tile hadronic calorimeter is separated into a large barrel and 
two smaller extended barrel cylinders, one on either side of 
the central baiTel. In the end-caps, copper/LAr technology is 
used for the hadronic end-cap calorimeters (HEC), covering 
the range 1.5 < |tj| < 3.2. The copper-tungsten/LAr forward 
calorimeters (FCal) provide both electromagnetic and hadronic 
energy measurements, extending the coverage to ji] | = 4.9. 

The trigger system consists of a hardware-based Level 1 
(LI) and a two-tier, software-based High Level Trigger (HLT). 
The LI jet trigger uses a sliding window algorithm with coarse- 
granularity calorimeter towers. This is then refined using jets 
reconstructed from calorimeter cells in the HLT. 



3 Monte Carlo simulation 

3.1 Event generators 

Data are compared to Monte Carlo (MC) simulations of 
jets with large transverse momentum produced via strong in- 
teractions described by Quantum Chromodynamics (QCD) in 
proton-proton collisions at a centre-of-mass energy of ^/s = 
7 TeV. The jet energy resolution is derived from several sim- 
ulation models in order to study its dependence on the event 
generator, on the parton showering and hadronisation models, 
and on tunes of other soft model parameters, such as those of 
the undelying event. The event generators used for this analysis 
are described below. 

1. Pythia 6.4 MClO tune: The event generator Pythia [15] 
simulates non-diffractive proton-proton collisions using a 
2 — > 2 matrix element at the leading order (LO) of the 
strong coupling constant to model the hard sub-process, 
and uses px-ordered parton showers to model additional ra- 
diation in the leading-logarithm approximation [16]. Mul- 
tiple parton interactions [17], as well as fragmentation and 
hadronization based on the Lund string model [18] are also 
simulated. The parton distribution function (PDF) set used 
is the modified leading-order MRST LO* set [19]. The pa- 
rameters used to describe multiple parton interactions are 
denoted as the ATLAS MClO tune [20]. This generator and 
tune are chosen as the baseline for the jet energy resolution 
studies. 

2. The Pythia Perugia2010 tune is an independent tune 
of Pythia to hadron collider data with increased final- 
state radiation to better reproduce the jet and hadronic event 
shapes observed in LEP and Tevatron data [21]. Parameters 
sensitive to the production of particles with strangeness and 
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related to jet fragmentation have also been adjusted. It is the 
tune favoured by ATLAS jet shape measurements [22]. 

3. The Pythia PARP90 modification is an independent sys- 
tematic variation of Pythia. The variation has been car- 
ried out by changing the parameter that controls the energy 
dependence of the cut-off, deciding whether the events are 
generated with the matrix element and parton-shower ap- 
proach, or the soft underlying event [23]. 

4. PythiaS [24] is based on the event generator PYTHIA and 
contains several modelling improvements, such as fully in- 
terleaved ;7t -ordered evolution of multiparton interactions 
and initial- and final-state radiation, and a richer mix of 
underlying-event processes. Once fully tested and tuned, it 
is expected to offer a complete replacement for version 6.4. 

5. The Herwigh-h- generator [25-28] uses a leading order 
2 — > 2 matrix element with angular-ordered parton showers 
in the leading-logarithm approximation. Hadronization is 
performed in the cluster model [29]. The underlying event 
and soft inclusive interactions use hard and soft multiple 
partonic interaction models [30]. The MRST LO* PDFs 
[19] are used. 

6. Alpgen is a tree-level matrix element generator for hard 
multi -parton processes (2 ri) in hadronic collisions [31]. 
It is interfaced to Herwig to produce parton showers in 
leading-logarithm approximation, which are matched to the 
matrix element partons with the MLM matching scheme 
[32]. Herwig is used for hadronization and Jimmy [33] 
is used to model soft multiple parton interactions. The LO 
CTEQ6L1 PDFs [34] are used. 



3.2 Simulation of the ATLAS detector 

Detector simulation is performed with the ATLAS simula- 
tion framework [35] based on Geant4 [36], which includes 
a detailed description of the geometry and the material of 
the detector. The set of processes that describe hadronic in- 
teractions in the Geant4 detector simulation are outlined in 
Refs. [37, 38]. The energy deposited by particles in the active 
detector material is converted into detector signals to mimic the 
detector read-out. Finally, the Monte Carlo generated events are 
processed through the trigger simulation of the experiment and 
are reconstructed and analysed with the same software that is 
used for data. 



3.3 Simulated pile-up samples 

The nominal MC simulation does not include additional 
proton-proton interactions (pile-up). In order to study its ef- 
fect on the jet energy resolution, two additional MC samples 
are used. The first one simulates additional proton-proton in- 
teractions in the same bunch crossing (in-time pile-up) while 
the second sample in addition simulates effects on calorimeter 
cell energies from close-by bunches (out-of-time pile-up). The 
average number of interactions per event is 1.7 (1.9) for the 
in-time (in-time plus out-of-time) pile-up samples, which is a 
good representation of the 2010 data. 
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4 Event and jet selection 



The status of each sub-detector and trigger, as well as re- 
constructed physics objects in ATLAS is continuously assessed 
by inspection of a standard set of distributions, and data-quality 
flags are recorded in a database for each luminosity block (of 
about two minutes of data-taking). This analysis selects events 
satisfying data-quality criteria for the Inner Detector and the 
calorimeters, and for track, jet, and missing transverse energy 
reconstruction [39]. 

For each event, the reconstructed primary vertex position is 
required to be consistent with the beamspot, both transversely 
and longitudinally, and to be reconstructed from at least five 
tracks with transverse momentum pj^^^ > 150 MeV associ- 
ated with it. The primary vertex is defined as the one with the 
highest associated sum of squared track transverse momenta 
L{pj^'^'^)^, where the sum runs over all tracks used in the ver- 
tex fit. Events are selected by requiring a specific OR combi- 
nation of inclusive single-jet and dijet calorimeter-based trig- 
gers [40,41]. The combinations are chosen such that the trigger 
efficiency for each pi bin is greater than 99%. For the lowest 
Pi bin (30^0 GeV), this requirement is relaxed, allowing the 
lowest-threshold calorimeter inclusive single-jet trigger to be 
used with an efficiency above 95%. 

Jets are reconstructed with the anti-^, jet algorithm [42] us- 
ing the FastJet software [43] with radius parameters R = 0.4 or 
R = 0.6, a four-momentum recombination scheme, and three- 
dimensional calorimeter topological clusters [44] as inputs. 
Topological clusters are built from calorimeter cells with a sig- 
nal at least four times higher than the root-mean-square (RMS) 
of the noise distribution (seed cells). Cells neighbouring the 
seed which have a signal to RMS -noise ratio > 2 are then iter- 
atively added. Finally, all nearest neighbour cells are added to 
the cluster without any threshold. 

Jets from non-collision backgrounds (e.g. beam-gas events) 
and instrumental noise are removed using the selection criteria 
outlined in Ref. [39]. 

Jets are categorized according to their reconstructed rapid- 
ity in four different regions to account for the differently instru- 
mented parts of the calorimeter: 

- Central region (IjI < 0.8). 

- Extended Tile Barrel (0.8 < \y\ < 1.2). 

- Transition region (1.2< \y\ <2.1). 

- End-Cap region (2.1 < \y\ < 2.8). 

Events are selected only if the transverse momenta of the two 
leading jets are above a jet reconstruction threshold of 7 GeV at 
the electromagnetic scale (see Section 5) and within \y\ < 2.8, 
at least one of them being in the central region. The analysis is 
restricted to \y\ < 2.8 because of the limited number of jets at 
higher rapidities. 

Monte Carlo simulated "particle jets" are defined as those 
built using the same jet algorithm as described above, but using 
instead as inputs the stable particles from the event generator 
(with a lifetime longer than 10 ps), excluding muons and neu- 
trinos. 
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5 Jet energy calibration 5.2 The Local Cluster Weighting (LCW) calibration 



Calorimeter jets are reconstructed from calorimeter energy 
deposits measured at the electromagnetic scale (EM-scale), the 
baseline signal scale for the energy deposited by electromag- 
netic showers in the calorimeter. Their transverse momentum 
is referred to . For hadrons this leads to a jet energy 

measurement that is typically 15-55% lower than the true en- 
ergy, due mainly to the non-compensating nature of the ATLAS 
calorimeter [45]. The jet response is defined as the ratio of 
calorimeter jet pj and particle jet pj, reconstructed with the 
same algorithm, and matched in 77 — space (see Section 9). 

Fluctuations of the hadronic shower, in particular of its 
electromagnetic content, as well as energy losses in the dead 
material lead to a degraded resolution and jet energy measure- 
ment compared to particles interacting only electromagneti- 
cally. Several complementary jet calibration schemes with dif- 
ferent levels of complexity and different sensitivity to system- 
atic effects have been developed to understand the jet energy 
measurements. The jet calibration is performed by applying 
corrections derived from Monte Carlo simulations to restore 
the jet response to unity. This is referred to as determining the 
jet energy scale (JES). 

The analysis presented in this article aims to determine the 
jet energy resolution for jets reconstructed using various JES 
strategies. A simple calibration, referred to as the EMh-JES cal- 
ibration scheme, has been chosen for the 2010 data [39]. It al- 
lows a direct evaluation of the systematic uncertainties from 
single-hadron response measurements and is therefore suitable 
for first physics analyses. More sophisticated calibration tech- 
niques to improve the jet resolution and reduce partonic flavour 
response differences have also been developed. They are the 
Local Cluster Weighting (LCW), the Global Cell Weighting 
(GCW) and the Global Sequential (GS) methods [39]. In ad- 
dition to these calorimeter calibration schemes, a Track-Based 
Jet Correction (TBJC) has been derived to adjust the response 
and reduce fluctuations on a jet-by-jet basis without changing 
the average jet energy scale. These calibration techniques are 
briefly described below. 



5.1 The EM+JES calibration 

For the analysis of the first proton-proton collisions, a sim- 
ple Monte Carlo simulation-based correction is applied as the 
default to restore the hadronic energy scale on average. The 
EMh-JES calibration scheme applies corrections as a function 
of the jet transverse momentum and pseudorapidity to jets re- 
constructed at the electromagnetic scale. The main advantage 
of this approach is that it allows the most direct evaluation of 
the systematic uncertainties. The uncertainty on the absolute jet 
energy scale was determined to be less than ±2.5% in the cen- 
tral calorimeter region (\y\ < 0.8) and ±14% in the most for- 
ward region (3.2 < \y\ < 4.5) for jets with pj > 30 GeV [39]. 
These uncertainties were evaluated using test-beam results, sin- 
gle hadron response in situ measurements, comparison with 
jets built from tracks, pj balance in dijet and 74-jet events, es- 
timations of pile-up energy deposits, and detailed Monte Carlo 
comparisons. 



The LCW calibration scheme uses properties of clusters to 
calibrate them individually prior to jet finding and reconstruc- 
tion. The calibration weights are determined from Monte Carlo 
simulations of charged and neutral pions according to the clus- 
ter topology measured in the calorimeter. The cluster properties 
used are the energy density in the cells forming them, the frac- 
tion of their energy deposited in the different calorimeter lay- 
ers, the cluster isolation and its depth in the calorimeter. Cor- 
rections are applied to the cluster energy to account for the en- 
ergy deposited in the calorimeter but outside of clusters and en- 
ergy deposited in material before and in between the calorime- 
ters. Jets are formed from calibrated clusters, and a final correc- 
tion is applied to the jet energy to account for jet-level effects. 
The resulting jet energy calibration is denoted as LCWh-JES. 

5.3 The Global Cell Weighting (GCW) calibration 

The GCW calibration scheme attempts to compensate for 
the different calorimeter response to hadronic and electromag- 
netic energy deposits at cell level. The hadronic signal is char- 
acterized by low cell energy densities and, thus, a positive 
weight is applied. The weights, which depend on the cell en- 
ergy density and the calorimeter layer only, are determined by 
minimizing the jet resolution evaluated by comparing recon- 
structed and particle jets in Monte Carlo simulation. They cor- 
rect for several effects at once (calorimeter non-compensation, 
dead material, etc.). A jet-level correction is applied to jets re- 
constructed from weighted cells to account for global effects. 
The resulting jet energy caUbration is denoted as GCWh-JES. 

5.4 The Global Sequential (GS) calibration 

The GS calibration scheme uses the longitudinal and trans- 
verse structure of the jet calorimeter shower to compensate for 
fluctuations in the jet energy measurement. In this scheme the 
jet energy response is first calibrated with the EMh-JES cali- 
bration. Subsequently, the jet properties are used to exploit the 
topology of the energy deposits in the calorimeter to charac- 
terize fluctuations in the hadronic shower development. These 
corrections are applied such that the mean jet energy is left un- 
changed, and each correction is applied sequentially. This cali- 
bration is designed to improve the jet energy resolution without 
changing the average jet energy scale. 

5.5 Track-based correction to the jet calibration 

Regardless of the inputs, algorithms and calibration meth- 
ods chosen for calorimeter jets, more information on the jet 
topology can be obtained from reconstructed tracks associated 
to the jet. Calibrated jets have an average energy response close 
to unity. However, the energy of an individual jet can be over- 
or underestimated depending on several factors, for example: 
the ratio of the electromagnetic and hadronic components of 
the jet; the fraction of energy lost in dead material, in either the 
inner detector, the solenoid, the cryostat before the LAr, or the 
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5 



cryostat between the LAr and the TileCal. The reconstructed 
tracks associated to the jet are sensitive to some of these ef- 
fects and therefore can be used to correct the calibration on a 
jet-by -jet basis. 

In the method referred to as Track-Based Jet Correction 
(TBJC) [45], the response is adjusted depending on the num- 
ber of tracks associated with the jet. The jet energy response 
is observed to decrease with jet track multiplicity mainly be- 
cause the ratio of the electromagnetic to the hadronic compo- 
nent decreases on average as the number of tracks increases. 
In effect, a low charged-track multiplicity typically indicates a 
predominance of neutral hadrons, in particular tt^s which yield 
electromagnetic deposits in the calorimeter with /? ~ 1 . A large 
number of charged particles, on the contrary, signals a more 
dominant hadronic component, with a lower response due to 
the non-compensating nature of the calorimeter {h/e < 1). The 
TBJC method is designed to be applied as an option in addi- 
tion to any JES calibration scheme, since it does not change 
the overall response, to reduce the jet-to-jet energy fluctuations 
and improve the resolution. 



where a{pj^i) — C7(/?t,2) — <y{pT), since both jets are in the 
same y region. 

If one of the two leading jets (j) is in the rapidity bin being 
probed and the other one (/) in a reference y region, it can be 
shown that the fractional jet pj resolution is given by 

^ I (,) = v'4ct2(A(,,,))-2ct2(A(,)) , (3) 

where ^) is measured in a topology with the two jets in dif- 
ferent rapidity regions and where (/) = (/,/) denotes both jets 
in the same y region. 

The back-to-back requirement is approximated by an az- 
imuthal angle cut between the leading jets, A(p{ji,j2) > 2.8, 
and a veto on the third jet momentum, p™^^'^^''^ < lOGeV, 
with no rapidity restriction. The resulting asymmetry distri- 
bution is shown in Fig. 1 for a pj = (px j + pj2)/2 bin of 
60GeV <pj< 80GeV, in the central region (\y\ < 0.8). Rea- 
sonable agreement in the bulk is observed between data and 
Monte Carlo simulation. 



6 In situ jet resolution measurement using 
the dijet balance method 

Two methods are used in dijet events to measure in situ the 
fractional jet pj resolution, (t(/?t) / Pt, which at fixed rapidity 
is equivalent to the fractional jet energy resolution, o{E)/E. 
The first method, presented in this section, relies on the ap- 
proximate scalar balance between the transverse momenta of 
the two leading jets and measures the sensitivity of this balance 
to the presence of extra jets directly from data. The second one, 
presented in the next section, uses the projection of the vector 
sum of the leading jets' transverse momenta on the coordinate 
system bisector of the azimuthal angle between the transverse 
momentum vectors of the two jets. It takes advantage of the 
very different sensitivities of each of these projections to the 
underlying physics of the dijet system and to the jet energy res- 
olution. 



6.1 Measurement of resolution from asymmetry 

The dijet balance method for the determination of the jet 
Pi resolution is based on momentum conservation in the trans- 
verse plane. The asymmetry between the transverse momenta 
of the two leading jets A(/5t,i,/?t,2) is defined as 



A(/?T,l,/'T,2) 



PT,1 —PT,2 
PTA + PT,2 



(1) 



where pj i and pj 2 refer to the randomly ordered transverse 
momenta of the two leading jets. The width (7(A) of a Gaus- 
sian fit to A{pj I, Pj 2) is used to characterize the asymmetry 
distribution and determine the jet pj resolutions. 

For events with exactly two particle jets that satisfy the hy- 
pothesis of momentum balance in the transverse plane, and re- 
quiring both jets to be in the same rapidity region, the relation 
between (t(A) and the fractional jet resolution is given by 



(T(A) 



y/a^{pT,i) + (y^{PT,2) 1 (y{PT) 



"D|T3 



— +- Data 2010 is = 7TeV 
I I Monte Carlo (PYTHIA) 



ATLAS 

Anti-k, R = 0.6 jets 
EM+JES calibration 
60 < < 80 GeV 

|y| < 0.8 




(PT,1 +PT,2) 



V2 Pt 



(2) 



Fig. 1: Asymmetry distribution as defined in Equation (1) for 
= 60 — 80 GeV and \y\ < 0.8. Data (points with error bars) 
and Monte Carlo simulation (histogram with shaded bands) are 
overlaid, together with a Gaussian fit to the data. The lower 
panel shows the ratio between data and MC simulation. The 
errors shown are only statistical. 



6.2 Soft radiation correction 

Although requirements on the azimuthal angle between the 
leading jets and on the third jet transverse momentum are de- 
signed to enrich the purity of the back-to-back jet sample, it is 
important to account for the presence of additional soft particle 
jets not detected in the calorimeter 

In order to estimate the value of the asymmetry for a pure 
particle dijet event, (y{pi)/pi = ^/2a{A) is recomputed al- 
lowing for the presence of an additional third jet in the sample 
for a series of cut-off threshold values up to 20 GeV. 



The cut on the third jet is placed at the EM-scale to be inde- 
pendent of caHbration effects and to have a stable reference for 
all calibration schemes. For each pj bin, the jet energy resolu- 



tions obtained with the different 



Pt3 



cuts are fitted with 
in order to 

estimate the expected resolution for an ideal dijet topology 



a straight line and extrapolated to ^'^^''^ — > 0, 
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of the form KsoftipT) = a + b / {\og p-r)^ , which was found to 
describe the distribution well, within uncertainties. The dif- 
ferences in the resolution due to other parameterizations were 
studied and treated as a systematic uncertainty, resulting in a 
relative uncertainty of about 6% (see Section 10). 



EM— scale 



a{pj) 



PT 



EM— scale 
PT3 



6.3 Particle balance correction 



The dependence of the jet pj resolution on the presence of a 
third jet is illustrated in Fig. 2. The linear fits and their extrapo- 
lations for a Pj bin of 60 < pj < 80 GeV are shown. Note that 
the resolutions become systematically broader as the pj^^^^^^'^ 
cut increases. This is a clear indication that the jet resolution 
determined from two-jet topologies depends on the presence of 
additional radiation and on the underlying event. 



The Pj difference between the two calorimeter jets is not 
solely due to resolution effects, but also to the balance between 
the respective particle jets. 
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Fig. 2: Fractional jet pj resolutions, from Equation 2, mea- 
sured in events with 60 < pj < 80 GeV and with third jet 
with PT less than 77™-''=^'% as a function of for data 

(squares) and Monte Carlo simulation (circles). The solid lines 
correspond to linear fits while the dashed lines show the ex- 
trapolations to p™^'*'^"''' = 0. The lower panel shows the ratio 
between data and MC simulation. The errors shown are only 
statistical. 

A soft radiation (SR) correction factor, Ksoftipi), is ob- 
tained from the ratio of the values of the linear fit at GeV 
and at 10 GeV: 



KsoftiPJ 









PT 


„EM— scale 

Pt.3 — 


^0 GeV 










„EM-scale 

Pt.3 - 


10 GeV 



(4) 



This multiplicative correction is applied to the resolutions 
extracted from the dijet asymmetry for p™"^'^^'*' < 10 GeV 
events. The correction varies from 25% for events with pj of 
50 GeV down to 5% for pj of 400 GeV. In order to limit the 
statistical fluctuations, KsoftipT) is fit with a parameterization 



The measured difference (left side) is decomposed into res- 
olution fluctuations (the first two terms on the right side) plus 
a particle-level balance (PB) term that originates from out-of- 
jet showering in the particle jets and from soft QCD effects. In 
order to correct for this contribution, the particle-level balance 
is estimated using the same technique (asymmetry plus soft ra- 
diation correction) as for calorimeter jets. The contribution of 
the dijet PB after the SR correction is subtracted in quadra- 
ture from the in situ resolution for both data and Monte Carlo 
simulation. The result of this procedure is shown for simulated 
events in the central region in Fig. 3. The relative size of the 
particle-level balance correction with respect to the measured 
resolutions varies between 2% and 10%. 
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Fig. 3: Fractional jet resolution obtained in simulation using 
the dijet balance method, shown as a function of pj, both be- 
fore (circles) and after the particle-balance (PB) correction (tri- 
angles). Also shown is the dijet PB correction itself (squares) 
and, in the lower panel, its relative size with respect to the frac- 
tional jet resolution. The errors shown are only statistical. 
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7 In situ jet resolution measurement using 
the bisector method 

7.1 Bisector rationale 

The bisector method [14] is based on a transverse balance 
vector, Py, defined as the vector sum of the momenta of the two 
leading jets in dijet events. This vector is projected along an 
orthogonal coordinate system in the transverse plane, (\//, Tj), 
where rj is chosen in the direction that bisects the angle formed 
by pj I and pj 2, '^(pn = ^ ^2- This is illustrated in Fig. 4. 




< 1 <■ 



The resolution is thus expressed in terms of calorimeter observ- 
ables only. The contribution from soft radiation and the under- 
lying event is minimised by subtracting in quadrature from 

If one of the leading jets (j) belongs to the rapidity region 
being probed, and the other one (/) to a previously measured 
reference y region, then 



Pi 



t2 calo 



-CT, 



2 calo 



U) 



/9^(|COS4012| 



a^{pi 



(6) 



PT,1 



The dispersions (7^ and Grj are extracted from Gaussian fits 
to the F^jf, and Pt.t? distributions in bins of px- There is no Aij) 
cut imposed between the leading jets, but it is implicitly limited 
by a p™^'*'^^''^ < 10 GeV requirement on the third jet, as dis- 
cussed in the next section. Figure 5 compares the distributions 
of Pj and Pj n between data and Monte Carlo simulation in 
the momentum bin 60 < pj < 80 GeV. The distributions agree 
within statistical fluctuations. The resolutions obtained from 
the Pj Y and Pj components of the balance vector are sum- 
marised in the central region as a function of pj in Fig. 6. As 
expected, the resolution on the Tj component does not vary with 
the jet Pj, while the resolution on the component degrades 
as the jet pj increases. 



7.2 Validation of the soft radiation isotropy with data 



Fig. 4: Variables used in the bisector method. The Tj-axis cor- 
responds to the azimuthal angular bisector of the dijet system 
in the plane transverse to the beam, while the i/A-axis is defined 
as the one orthogonal to the Tj-axis. 

For a perfectly balanced dijet event, Pj = 0. There are of 
course a number of sources that give rise to significant fluc- 
tuations around this value, and thus to a non-zero variance of 
its ij/ and tj components, denoted (7^ and a^, respectively. At 

particle level, Pj^^^ receives contributions mostly from initial- 
state radiation. This effect is expected to be isotropic in the 
(i/, T]) plane, leading to similar fluctuations in both compo- 
nents, CT^"'^' = ajj"'^. 

The validity of this assumption, which is at the root of the 
bisector method, can be checked with Monte Carlo simulations 
and with data. The precision with which it can be assessed 
is considered as a systematic uncertainty (see Section 7.2). 
The 1/ component has greater sensitivity to the energy reso- 
lution because Pr,i/ is the difference between two large trans- 
verse momentum components while Pj ^ is the sum of two 
small components. Effects such as contamination from 3 -jet 
events or final-state radiation not absorbed in the leading jets 
by the clustering algorithm could give rise to a a^"' > ct}^"'. At 

^2 calo 



Figure 7 shows the width of the \j/ and tj components of Pj- 
as a function of the p™^^'^^''^ cut, for anti-A:, jets with R — 0.6. 
The two leading jets are required to be in the same rapidity re- 
gion, \y\ < 0.8, while there is no rapidity restriction for the third 
jet. As expected, both components increase due to the contri- 
bution from soft radiation as the pj $ cut is increased. Also 
shown as a function of the p™^'*'^^''^ cut is the square-root of 
the difference between their variances, which yields the frac- 
tional momentum resolution when divided by 2 (pj) (cos A^). 

It is observed that the increase of the soft radiation contri- 
bution to CT^^'° and CT^^'° cancels in the squared difference and 
that it remains almost constant, within statistical uncertainties, 
up to pf^-'"''^^ ~ 20 GeV for pj between 160-260 GeV. The 
same behaviour is observed for other pj ranges. This cancella- 
tion demonstrates that the isotropy assumption used for the bi- 
sector method is valid over a wide range of choices of p™^^'^^'*' 
without the need for requiring an explicit A (j) cut between the 
leading jets. The precision with which it can be ascertained in 

part 



(7n is taken conservatively as a systematic 



situ that (7^^" 

uncertainty on the method, of about 4 — 5% at 50 GeV (see 
Section 10). 



calorimeter level, (7^ calo expected to be significantly larger 8 Performance for the EM+JES Calibration 



than (7^ mostly because of the jet energy resolution. 

If both jets belong to the same y region, such that they have 
the same average jet energy resolution, it can be shown that 



(t2 calo _(y2 calo 



Pt V2ptV(|cos4(/)i2|)' 



(5) 



The performances of the dijet balance and bisector meth- 
ods are compared for both data and Monte Carlo simulation as 
a function of jet pj for jets reconstructed in the central region 
with the anti-^, algorithm with and using the EMh-JES 

calibration scheme. The results are shown in Fig. 8. The resolu- 
tions obtained from the two independent in situ methods are in 
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Fig. 5: Distributions of the Pj (top) and fx.r] (bottom) compo- 
nents of the balance vector Py, for pj — 60 — 80 GeV. The data 
(points with error bars) and Monte Carlo simulation (histogram 
with shaded bands) are overlaid. The lower panel shows the ra- 
tio between data and MC simulation. The errors shown are only 
statistical. 

good agreement with each other within the statistical uncertain- 
ties. The agreement between data and Monte Carlo simulation 
is also good with some deviations observed at low pj. 

The resolutions for the three jet rapidity bins with \y\ > 0.8, 
the Extended Tile Barrel, the Transition and the End-Cap re- 
gions, are measured using Eqs. 3 and 6, taking the central re- 
gion as the reference. The results for the bisector method are 
shown in Fig. 9. Within statistical errors the resolutions ob- 
tained for data and Monte Carlo simulation are in agreement 
within ±10% over most of the px-range in the various regions. 

Figure 9 shows that dependences are well described by fits 
to the standard functional form expected for calorimeter-based 
resolutions, with three independent contributions, the effective 
noise (A^), stochastic (S) and constant (C) terms. 



Pj 



N 

PT 



c. 



Pi 



(7) 



The term is due to external noise contributions that are not (or 
only weakly) dependent on the jet pj, and include the electron- 



> 
<3 



50 



40 



30- 



20 



10 




Bisector Method 
|y|<0.8 



Anti-k, R= 0.6 jets 
EM+JES calibration 



Pt Data 
P * Data 

P^jl Monte Carlo (PYTHIA) 
PjI Monte Carlo (PYTHIA) - 



-A-- 



-H- 



ATLAS 

H \ h 





-•a H js 

Data 2010 \'s = 



: 7TeV 



:j:.::J:. 



30 40 50 



100 



200 300 
(p^ +p^^)/2 (GeV) 

Fig. 6: Standard deviations of Pj ^ and Pt.tj. the components of 
the balance vector, as a function of pj. The lower panel shows 
the ratio between data and MC simulation. The errors shown 
are only statistical. 
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Fig. 7: Standard deviations cr^"'", (7^^^° and [(cr^ - (7^) 
as a function of the upper p™"^'^^'*' cut, for R = 0.6 anti-fc, jets 
with Pj = 160 — 260 GeV. The errors shown are only statistical. 



ics and detector noise, and contributions from pile-up. It is ex- 
pected to be significant in the low-pT region, below ^30 GeV. 
The C term encompasses the fluctuations that are a constant 
fraction of the jet pj, assumed at this early stage of data-taking 
to be due to real signal lost in passive material (e.g. cryostats 
and solenoid coil), to non-uniformities of response across the 
calorimeter, etc. It is expected to dominate the high-y^x region, 
above 400 GeV. For intermediate values of the jet pj, the sta- 
tistical fluctuations, represented by the S term, become the lim- 
iting factor in the resolution. With the present data sample that 
covers a restricted pj range, 30 GeV < pj < 500 GeV, there is 
a high degree of correlation between the fitted parameters and it 
is not possible to unequivocally disentangle their contributions. 



ATLAS collaboration: Jet Energy Resolution with the ATLAS detector 



9 



H 0.25 

Q. 
Cl 

■o 0.2 
0.15 
0.1 
0.05 



5? 


20 




O 










O 







lb 








■H 
Q 




-20 





>!irz.>!is 



Data 2010 Vs = 7TeV 



Anti-k, R = 0.6 jets 
EM+JES calibration 

|y,l<0.8 







Bisector: Data 

Dijet Balance: Data 

Bisector: Monte Carlo (PYTHIA) 

Dijet Balance: Monte Carlo (PYTHIA) 



:4:::::a:::::::^:: 



t % 



30 40 50 60 70 80 100 



200 300 400 500 

(P,,+P,,)/2 (GeV) 

Fig. 8: Fractional jet pj resolution for the dijet balance and bi- 
sector methods as a function of pj. The lower panel shows the 
relative difference between data and Monte Carlo results. The 
dotted lines indicate a relative difference of ±10%. Both meth- 
ods are found to be in agreement within 10% between data and 
Monte Carlo simulation. The errors shown are only statistical. 

9 Closure test using Monte Carlo simulation 



10 Jet energy resolution uncertainties 
10.1 Experimental uncertainties 

The squares (circles) in Fig. 11 show the experimental 
relative systematic uncertainty in the dijet balance (bisector) 
method as a function of pj. The different contributions are dis- 
cussed below. The shaded area corresponds to the larger of the 
two systematic uncertainties for each pj bin. 

For the dijet balance method, systematic uncertainties take 
into account the variation in resolution when applying different 
40 cuts (varied from 2.6 to 3.0), resulting in a 2-3% effect for 
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The Monte Carlo simulation expected resolution is de- 
rived considering matched particle and calorimeter jets in the 
event, with no back-to-back geometry requirements. Match- 
ing is done in Tj - space, and jets are associated if 
AR = \J (4tj)2 + (40)2 < 0.3. The jet response is defined as ^ o.3 

in bins of /^ip"", where /9Sj?'° and pp^^^ correspond 
to the transverse momentum of the reconstructed jet and its 
matched particle jet, respectively. The jet response distribution 
is modelled with a Gaussian fit, and its standard deviation is 
defined as the truth jet p-\ resolution. 

The Monte Carlo simulation truth jet p-\ resolution is com- 
pared to the results obtained from the dijet balance and the 
bisector in situ methods (applied to Monte Carlo simulation) 
in Fig. 10. The agreement between the three sets of points is 
within 10%. This result confirms the validity of the physical as- 
sumptions discussed in Sections 6 and 7 and the inference that 
the observables derived for the in situ MC dijet balance and 
bisector methods provide reliable estimates of the jet energy 
resolution. The systematic uncertainties on these estimates are 
of the order of 10% (15%) for jets with R = 0.6 (R = 0.4), and 
are discussed in Section 10. 
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Fig. 9: Fractional jet pi resolution as a function of pj for anti- 
kt with R = 0.6 jets in the Extended Tile Barrel (top). Transi- 
tion (center) and End-Cap (bottom) regions using the bisector 
method. In the lower panel of each figure, the relative differ- 
ence between the data and the MC simulation results is shown. 
The dotted lines indicate a relative difference of ±10%. The 
errors shown are only statistical. 



Pj = 30-60 GeV, and when varying the soft radiation correc- 
tion modelling, which contributes up to 6% at pj « 30 GeV. 
For the bisector method, the relative systematic uncertainty is 
about 4-5%, and is derived from the precision with which the 



assumption that cr^^ 
can be verified 



part 



when varying the pj 



EM— scale 



cut 



The contribution from the JES uncertainties [39] is 1-2%, 
determined by re-calculating the jet resolutions after varying 
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Fig. 10: Comparison between the Monte Carlo simulation truth 
jet pT resolution and the results obtained from the bisector and 
dijet balance in situ methods (applied to Monte Carlo simu- 
lation) for the EM+JES calibration, as a function of pj. The 
lower panel of the figure shows the relative difference, obtained 
from the fits, between the in situ methods and Monte Carlo 
truth results. The dotted lines indicate a relative difference of 
±10%. The errors shown are only statistical. 

the JES within its uncertainty in a fully correlated way. The 
resolution has also been studied in simulated events with added 
pile-up events (i.e. additional interactions as explained in Sec- 
tion 3.3), as compared to events with one hard interaction only. 
The sensitivity of the resolution to pile-up is found to be less 
than 1% for an average number of vertices per event of 1.9. 

In summary, the overall relative uncertainty from the in situ 
methods decreases from about 7% at pj =30 GeV down to 
4% at Pj — 500 GeV. Figure 1 1 also shows in dashed lines 
the absolute value of the relative difference between the two in 
situ methods, for both data and Monte Carlo simulation. They 
are found to be in agreement within 4% up to 500 GeV, and 
consistent with these systematic uncertainties. 
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Fig. 11: The experimental systematic uncertainty on the dijet 
balance (squares) and bisector (circles) methods as a function 
of Pj, for jets with \y\ < 0.8. The absolute value of the rela- 
tive difference between the two methods in each pj bin is also 
shown for data and for Monte Carlo simulation (dashed lines). 



10.2 Uncertainties due to the event modelling in the 
Monte Carlo generators 

The expected jet pj resolution is calculated for other Monte 
Carlo simulations in order to assess its dependence on different 
generator models (Alpgen and HERWIG++), PYTHIA tunes 
(PERUGIA2010), and other systematic variations (PARP90; see 
Sec. 3.1). Differences between the nominal Monte Carlo simu- 
lation and PythiaS [24] have also been considered. These ef- 
fects, displayed in Fig. 12, never exceed 4%. Although they are 
not relevant for the in situ measurements of the jet energy res- 
olution themselves, physics analyses sensitive to the expected 
resolution have to consider a systematic uncertainty from event 
modelling estimated from the sum in quadrature of the differ- 
ent cases considered here. This is shown by the shaded area in 
Fig. 12 and found to be at most 5%. 
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Fig. 12: Systematic uncertainty due to event modelling in 
Monte Carlo generators on the expected jet energy resolution as 
a function of pj, for jets with \y\ < 0.8. The reference is taken 
from Pythia MClO and other event generators are shown 
as solid triangles (HERWIG++) and open circles (Alpgen). 
Solid squares (PYTHIA PERUGIA2010) and inverted triangles 
(Pythia PARP90) summarize differences coming from differ- 
ent tunes and cut-off parameters, respectively. Open squares 
compare the nominal simulation with PythiaS. 



10.3 Uncertainties on the measured resolutions 

The uncertainties in the measured resolutions are dominated 
by the systematic uncertainties, which are shown in Table 1 
as a percentage of the resolution for the four rapidity regions 
and the two jet sizes considered, and for characteristic ranges, 
low (- 50 GeV), medium (- 150 GeV) and high (- 400 GeV) 
Pj. The results are similar for the four calibration schemes. 
The dominant sources of systematic uncertainty are the closure 
and the data/MC agreement. The closure uncertainty (see Sec- 
tion 9), defined as the precision with which in simulation the 
resolution determined using the in situ method reproduces the 
truth jet resolution, is larger for R ~ 0.4 than for R — 0.6, de- 
creases with Pj, and is basically independent of the rapidity. 
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Jet 


Rapidity 


Total systematic uncertainty 


radius 


range 


1 J^\I7 nn-i 

J_,UW jJ^ 


iViCQ yyy 


mgn p-Y 




n <' lul <^ n s 
^ — l/l ^ U.o 


12% 


10% 


11% 


R = 0.6 


0.8 < \y\ < 1.2 


12% 


10% 


13% 




1.2 < \y\ <2.1 


14% 


12% 


14% 




2.1 < \y\ <2.8 


15% 


13% 


18% 




o< bl <0.8 


17% 


15% 


11% 


R = OA 


0.8 < \y\ < 1.2 


20% 


18% 


14% 




1.2 < \y\ <2.1 


20% 


18% 


14% 




2.1 < \y\ <2.8 


20% 


18% 


18% 



Table 1 : Relative systematic uncertainties at low 50 GeV), 
medium (- 150 GeV) and high (- 400 GeV) pj, for the four 
rapidity regions and the two jet radii studied. The uncertainties 
are similar for the four calibration schemes. 



The data/MC agreement uncertainty is observed to be indepen- 
dent of R, larger at low and high pj than at medium pj, and to 
grow with rapidity because of the increasingly limited statisti- 
cal accuracy with which checks can be performed to assess it. 
Other systematic uncertainties are significantly smaller. They 
include the validity of the soft radiation hypothesis, the jet en- 
ergy scale uncertainty and the dependence on the number of 
pile-up interactions. The uncertainty due to event modelling is 



not included, as it does not contribute to the in situ measure- 
ment itself. 

The systematic uncertainties in Table 1 for jets with R = 0.4 
are dominated by the contribution from the closure test. They 
decrease with pj and are constant for the highest three rapidity 
bins. They are also consistently larger than for the /? = 0.6 case. 
The systematic uncertainties for jets with R = 0.6 receive com- 
parable contributions from closure and data/MC agreement. 
They tend to increase with rapidity and are slightly lower in 
the medium pj range. The uncertainty increases at high pj for 
the end-cap, 2.1 < \y\ < 2.8, because of the limited number of 
events in this region. 



11 Jet energy resolution for other 
calibration schemes 

The resolution performance for anti-fc, jets with R = 0.6 
reconstructed from calorimeter topological clusters for the Lo- 
cal Cluster Weighting (LCWh-JES), the Global Cell Weight- 
ing (GCWh-JES) and the Global Sequential (GS) calibration 
strategies (using the bisector method) is presented in Fig. 13 
for the Central, Extended Tile Barrel, Transition and End-Cap 
regions. The top part shows the resolutions determined from 
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Fig. 13: Fractional jet pj resolutions as a function of pj for anti-A:, jets with R ~ 0.6 with \y\ < 0.8 (top left), 0.8 < [y] < 1 .2 (top 
right), 1 .2 < \y\ < 2. 1 (bottom left) and 2.1 < \y\ < 2.8 (bottom right), using the bisector in situ method, for four jet calibration 
schemes: EMh-JES, Local Cluster Weighting (LCWh-JES), Global Cell Weighting (GCWh-JES) and Global Sequential (GS). The 
lower panels show the relative difference between data and Monte Carlo simulation results. The dotted lines indicate relative 
differences of ±10%. The errors shown are only statistical. 
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data, whereas the bottom part compares data and Monte Carlo 
simulation results. The relative improvement in resolution with 
respect to the EM+JES calibrated jets is comparable for the 
three more sophisticated calibration techniques. It ranges from 
10% at low pj up to 40% at high pj for all four rapidity re- 
gions. 

Figure 14 displays the resolutions for the two in situ meth- 
ods applied to data and Monte Carlo simulation for \y\ < 0.8 
(left plots). It can be observed that the results from the two 
methods agree, within uncertainties. The Monte Carlo simula- 



tion reproduces the data within 10%. The figures on the right 
show the results of a study of the closure for each case, where 
the truth resolution is compared to that obtained from the in situ 
methods applied to Monte Carlo simulation data. The agree- 
ment is within 10%. Overall, comparable agreement in reso- 
lution is observed in data and Monte Carlo simulation for the 
EMh-JES, LCWh-JES, GCWh-JES and GS calibration schemes, 
with similar systematic uncertainties in the resolutions deter- 
mined using in situ methods. 
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Fig. 14: Fractional jet pj resolutions as a function of pj for anti-A:; jets with R = 0.6 for the Local Cluster Weighting (LCWh-JES), 
Global Cell Weighting (GCWh-JES) and Global Sequential (GS) calibrations. Left: Comparison of both in situ methods on data 
and MC simulation for \y\ < 0.8. The lower panels show the relative difference. Right: Comparison between the Monte Carlo 
simulation truth jet pj resolution and the final results obtained from the bisector and dijet balance in situ methods (applied to 
Monte Carlo simulation). The lower panels show the relative differences, obtained from the fits, between the in situ methods and 
Monte Carlo truth results. The dotted lines indicate relative differences of ±10%. The errors shown are only statistical. 
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12 Improvement in jet energy resolution 
using tracks 

The addition of tracking information to the calorimeter- 
based energy measurement is expected to compensate for the 
jet-by-jet fluctuations and improve the jet energy resolution 
(see Section 5.5). The performance of the Track-Based Jet Cor- 
rection method (TBJC) is studied by applying it to both the 
EMh-JES and LCWh- JES calibration schemes, in the central 
region. The measured resolution for anti-^, jets with R = 0.6 
(R — 0.4) is presented as a function of the average jet trans- 
verse momentum in the top (bottom) plot of Fig. 15. 

The relative improvement in resolution due to the addition 
of tracking information is larger at low pj and more impor- 
tant for the EMh-JES calibration scheme. It ranges from 22% 
(10%) at low pj to 15% (5%) at high pj for the EMh-JES 
(LCW+JES) calibration. For pj < 70 GeV, jets calibrated with 
the EM+JES+TBJC scheme show a similar performance to 
those calibrated with the LCW+JES+TBJC scheme. Overall, 
jets with LCW+JES+TBJC show the best fractional energy res- 
olution over the full pj range. 
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Fig. 15: Top: Fractional jet pj resolutions as a function pj, 
measured in data for anti-fc, jets with R = 0.6 (top) and R = 
0.4 (bottom) and for four jet calibration schemes: EM+JES, 
EM+JES+TBJC, LCW+JES and LCW+JES+TBJC. The lower 
panel of the figure shows the relative improvement for the 
EM+JES+TBJC, LCW+JES and LCW+JES+TBJC calibra- 
tions with respect to the EM+JES jet calibration scheme, used 
as reference (dotted line). The errors shown are only statistical. 



13 Summary 



The jet energy resolution for various JES calibration schemes 
has been measured using two in situ methods with a data sam- 
ple corresponding to an integrated luminosity of 35 pb^' col- 
lected in 2010 by the ATLAS experiment at = 7 TeV. 

The Monte Carlo simulation describes the jet energy res- 
olution measured in data within 10% for jets with pj values 
between 30 GeV and 500 GeV in the rapidity range \y\ < 2.8. 

The resolutions obtained applying the in situ techniques to 
Monte Carlo simulation are in agreement within 10% with the 
resolutions determined by comparing jets at calorimeter and 
particle level. Overall, the results measured with the two in situ 
methods have been found to be consistent within systematic 
uncertainties. 
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